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Abstract 

This paper gives a note on an application of the enclosure method to an inverse 
obstacle scattering problem governed by the Helmholtz equation in two dimensions. 
It is shown that one can uniquely determine the convex hull of an unknown sound- 
hard polygonal obstacle from the trace of the total wave that was exerted by a 
single point source onto a known circle surrounding the obstacle provided the source 
is sufficiently far from the obstacle. The result contains a formula that extracts 
the value of the support function of the obstacle at a generic direction. Some 
other applications to thin obstacles, obstacles in a layered medium and the far-field 
equation in the linear sampling method are also included. 
AMS: 35R30 
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1 Introduction 

The enclosure method was introduced in [12] for inverse boundary value problems for ellip- 
tic equations which are motivated by the possibility of applications to electrical impedance 
tomography, diffraction tomography, etc.. Therein the observation data are formulated by 
using the Dirichlet-to-Neumann map (or Neumann-to-Dirichlet map) associated with the 
governing equation of a 'signal' propagating inside the medium. It aims at extracting in- 
formation about the location and shape of unknown discontinuity embedded in a known 
reference medium that gives an effect on the propagation of the signal, such as an obstacle, 
inclusion, crack, etc. from data observed on the boundary of the medium. Now we have 
many applications of this method, see, e.g., [11, 24, 17, 23, 25]. 

In [10] it was shown that, in a simplified situation a single set of the Dirichlet and 
Neumann data gives information about the convex hull of unknown discontinuity. It was 
the starting point of the single measurement version of the enclosure method and we have 
already many applications, e.g., [16, 13, 14, 18, 19, 15, 20, 21, 22]. 

This paper is closely related to [16]. Therein we considered an inverse obstacle scat- 
tering problem of acoustic wave in two dimensions. The problem is to reconstruct a 
two-dimensional obstacle from the Cauchy data on a circle surrounding the obstacle of 
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the total wave field generated by a single incident plane wave with a fixed wave number. 
Let us make a review of one of the results therein. 

We consider a polygonal obstacle denoted by D, that is, D C R 2 takes the form 
Di U • • ■ U D m with 1 < m < oo where each Dj is open and a polygon; IX,- D Dy = if 

./ / ./'• 

The total wave field u outside obstacle D takes the form u(x; d, k) = e tkx ' d + w(x) with 
k > 0, d E S 1 and satisfies 

Am + k 2 u = in R 2 \ D, 

^ = 0ondD, 
av 

(Ow \ 
— ikw = 0, 
or J 

where r = \x\ and v denotes the unit outward normal relative to D. The last condition 
above is called the Sommerfeld radiation condition. 

Let B R be an open disc with radius R centered at a fixed point satisfying D C B R . 
We assume that B R is known. Our data are u = u(-;d,k) and dujdv on dB R for a 
fixed d and k, where v is the unit outward normal relative to B R . Let u and oo 1 - be two 
unit vectors perpendicular to each other. We always choose the orientation of cu -1 and ui 
coincides with that of e 1 and e 2 and thus u 1 - is unique. 

We make use of the special complex exponential solution of the Helmholtz equation 
(A + k 2 )v = in R 2 : 

v T {x;u) = e *-(™W^Wu^ x £ r2) 

where r > is a parameter. 

Recall the support function of D: h^ius) = sup^^, x ■ u. We say that uj is regular with 
respect to D if the set dD H {x e~R 2 \ x ■ u = hoioj)} consists of only one point. 

Define 

J(t; w, tZ, fc) = jT^ _ ^7 M ) 

Theorem 1.1 ([16]). Assume that oj is regular with respect to D. Then the formula 

lim — log |/(r; a;, d, k)\ — }id(w), 

is valid. Moreover, we have the following: 

if t > ho(oj) , then lim e~ r< \I(t; oj, d, k)\ = 0; 

if t < h D {oj) , then lim e~ Tt |/(r; oj, d, k)\ = oo. 

In [19] a similar formula has been established by using the far-field pattern Frj(ip, d; k), 
(f G S 1 of scattered wave w = u — e lkx ' d for fixed d and k which determines the leading 
term of the asymptotic expansion of w as r — > oo in the following sense: 

w(rip) ~ —=F D (ip,d; k). 
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Moreover, therein instead of volumetric obstacle, similar formulae for thin sound-hard 
obstacle (or screen ) have also been established with two incident plane waves. 

In this section, we describe another inverse obstacle scattering problem in which a point 
source located within a finite distance from an unknown obstacle generates a scattered 
wave and one measures the total wave on a known circle surrounding an unknown obstacle. 
One can see this type of problem in, e.g., a mathematical formulation of microwave 
tomography [31], subsurface radar [5], etc.. 

Let y e R 2 \ D. Let E = E D (x,y) be the unique solution of the scattering problem: 

(A + k 2 )E = inR 2 \D, 
—E = -—%{■ ,y) ondD, 

(dE \ 
— ikE = 0, 
or J 



where 



*o(x,y) = l l H^\k\x-y\) 



and Hq 1 ^ denotes the Hankel function of the first kind [27]. 

The total wave outside D exerted by the point source located at y is given by the 
formula: 

<b D {x, y) = $ (x, y) + E D (x, y),xeR 2 \D. 

Inverse Problem. Let R\ > R. Fix k > and y e dB^. Extract information about 
the location and shape of D from <& D (x,y) given at all x G dB R . 

The aim of this paper is to show that the single measurement version of the enclosure 
method still works for this problem. 

Define 

J(r;u,y,k) = ^-§ D (x,y) -v t (x;uj) - -^-v t (x;uj) ■ $ D (a;,y)^ dS(x). 

The first result of this paper is as follows. 
Theorem 1.2. Assume that uj is regular with respect to D and that 

diamD < dist(D,dB Rl ). (1.1) 

It holds that ^ 

^lin^ - log | J(r;w, j/, k) \ = h D {u). 

Moreover, we have the following: 

if t > h D {oj) , then hm e~ r< | J(r; u, y, k)\ = 0; 

if t < hoioj) , then lim e~ r< | J(r; u, y, k)\ = oo. 



3 



It should be pointed out that (d/du)^o(x,y) for x G 8Br can be computed from 
$d(£, y) for x G dB R by solving the exterior Dirichlet problem for the Helmholtz equation: 



(A + A; 2 )£ = inR 2 \B fl , 
£ = $ D (-,y)-$ (-,y) ondfl fl , 

(<9-£/ ~ \ 
ikE I =0. 
c/r / 

The computation formula is 

(9 d d 

—$ D (x,y) = —$ (x,y) + q^E(x) ondB R . 

Needless to say, this kind of remak works also for dujdv in Theorem 1.1. 

Condition (1.1) can be satisfied if R\ is sufficiently large compared with R. It is not 
known whether condition (1.1) can be dropped completely. To suggest a possibility next 
we present a partial result which does not employ (1.1). 

For the description of the second result we introduce special scattered and total fields. 
Given d G S 1 choose d G S l in such a way that d L = d. Let x G dD and w = 
w(x] —d, k,x ) be the unique solution of the scattering problem: 

(A + k 2 )w = inR 2 \D, 

(fill (J 

_ = -_{(*„-.). tfe-*-} on M>, (12) 
lim \fr- l^r- — ikw] = 0. 

r — >oo \ Qf J 

Define 

u(x; -d, k, x ) — (x - x) ■ d e ~ tkx ' d + w(x; -d, k, x ). (1.3) 

Note that the function x i — > (xq — x) • <0 e - lkx - d satisfies the Helmholtz equation in the 
whole plane and the radiation condition for w(x; —d, k,x ) yields that 

u(x; -d, k, x ) = (x - x) ■ $ e - ikx - d + O^ 1 / 2 ) (1.4) 

as r — > oo. The second result of this paper is as follows. 

Theorem 1.3. Let oo be regular with respect to D and xq G 3D be the point satisfying 
x -oo = hoijjj)- Assume thaty G dB Rl satisfies u(y; —d±, k,Xo) ^ oru(y; —d 2 , k,x ) ^ 0, 
where di,d 2 G S 1 are directions of two sides of D that meet at xq. Then all the conclusions 
in Theorem 1.2 are valid. 

Note that in this theorem (1.1) is not assumed at the price of introducing another 
implicit restriction on the location of y relative to D. From (1.4) we see that given 
5 > and d G S* 1 if Rq > is sufficiently large, for all R± > Rq and y G dB Rl with 
\(xq — y) • i?| > 5 it holds that u(y; —d, k, x ) ^ with d = -i? -1 . We leave the problem of 
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removing the condition completely as an open problem. See also Remark 2.2 in Section 
2. 

Since the set of all oj which is not regular with respect to given D is finite and the 
support function of D is continuous, as a corollary of Theorem 1.2 we have uniqueness of 
determining the convex hull of D from $d(:e, y) given at all x G dB R for a single y G dB Rl 
and k > provided (1.1). It seems that this type of uniqueness with a single point source 
had not appeared in the previous study. See [7] and references therein for uniqueness of 
a polygonal obstacle with the far-field pattern of the scattered wave exerted by a single 
plane wave. It seems that their argument heavily depends on the fact that the total wave 
approaches the plane wave at infinity. This is not true for the total wave exerted by a 
point source since we have, as r — > oo, 

qItt/4 ikr 

$ D (rd,y)~ -j=— u ( y -,-d,k). (1.5) 

See [30, 9, 28] for the derivation for D with a smooth boundary. In our case dD is not 
smooth; however, a minor modification of the proof still works. Equation (1.5) means 
that the far-field pattern of $£>(^,y) as a function of x is given by u(y; —d, k) multiplied 
by a known constant. This formula has been used in the probe method [9] and the singular 
sources method [29]. 

Theorem 1.2 together with (1.5) yields 

Corollary 1.1. Let R x > R and D C B R . Assume that D satisfies (1.1). Fix k > and 
y G dB Rl . One can uniquely determine the convex hull of D from the data u(y; —d,k) 
given at all d G S 1 . 

Proof. The proof is divided into three steps. 

(i) Use (1.5) to compute the far-field pattern of $d(x, y) from u(y;—d,k) given at all 
d G S 1 . 

(ii) Use, e.g., the point source method [28] to compute $£>(x,y) together with its normal 
derivative for x G OBr from the far-field pattern of $d{x, y)- 

(iii) Use Theorem 1.2 to compute ho{oj) for a generic u from &£>(x,y) together with its 
normal derivative for x G 8Br. 

□ 

Summing up, we obtained two procedures for estimating the convex hull of an unknown 
sound-hard polygonal obstacle by using two types of the data. 
The first type of the data is given by the following process: 

(A) produce the total wave by a fixed point source located outside a known circle 
surrounding an unknown obstacle and observe the wave at all points on the circle. 

The second is as follows: 

(B) produce the total waves by incident plane waves for all directions and observe the 
waves at a fixed point outside a known circle surrounding an unknown obstacle. 

Note that these are different from the reciprocity principle ([4]) which is the identity 
Fo((p,d;k) = Fr)(—d,—(p;k) since in this identity the incident wave is always a plane 
wave and one observes the scattered wave at infinity. 

A brief outline of this paper is as follows. Theorems 1.2 and 1.3 are proved in Section 
2. Both proofs have a common starting point with the proof of Theorem 1.1 which we 
recall before describing subsections 2.1 and 2.2. In those subsections we complete the 
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proof of Theorems 1.2 and 1.3. In the last section three other applications are given. 
Two of them are concerned with some extensions to thin obstacles and obstacles in a 
layered medium. In the last of the applications we consider the far-field equation which 
plays the central role in the linear sampling method [3]. We show that a modification of 
the argument for the proof of Theorem 1.2 gives unsolvabilty of the far- field equation for 
polygonal obstacles. 



2 Proof of Theorems 1.2 and 1.3 

First we follow the argument for the proof of Theorem 1.1 (see also [18]). For simplicity 
of notation we set u(x) = Qd{x,v)- 

Let xo denote the single point of the set {x \ x ■ u> = ho(^)} H dD. x has to be a 
vertex of Dj for some j. In what follows we denote by B R (x ) the open disc with radius 
R centered at x . Let 6 denote the outside angle of D at i - @ satisfies n < O < 2tc 
since u is regular with respect to D. 

If one chooses a sufficiently small i] > 0, then one can write 

B 2v (x ) n (B Rl \ D) = {x + r(cos 9a + sin 9a L ) | < r < 2rj, < 9 < 0}, 

b v {x ) n dD = r p u r q u {x } 

where a = cos p co ± + sin p to, a L = — sin p to 1 - + cos p to; —tt < p < 0; T p = {x + ra \ < 
r < r]}, T q = {x + r(cos Qa + sin 6a 1 ) | < r < r/}. Note that the orientation of a, er 1 
coincides with that of e±, e 2 . See also Figure 1 of [10]. 

The quantity — p means the angle between two vectors oo 1 - and a. p satisfies > 
7r + (— p). Set q = — 2n + p. Then we have — n < q < p < and the expression 

T p = {xo + r(cos poj ± + sin pcu) \ < r < i]}, 

Tg = {^o + ^(cos gw 1 + sin quo) \ < r < r]}. 

This is the meaning of p and q. 
We set 

u(r, 9) = u(x), x = x + r(cos 9 a + sin 9 or 1 ). 
The u can be expanded as 

oo 

u(r, 9) = aiJ (kr) + OL n J\ n {kr) cos\ n 9, < r < 77, < 9 < 0, 

n=2 

where the A n describes the singularity of w as r — )> and in this case explicitly given by 
the formula A„ — (n — l)7r/0, J An stands for the Bessel function of order \ n . 

One of key points is introducing a new parameter s instead of r by the equation 
s = Vr 2 + k 2 + t, we obtain, as s — > 00, the complete asymptotic expansion 

J(r; A;) e"*-^^ ~ - i V 6 fc (2.1) 

n=2 6 
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where K n are constants given by the formula 

K n = e ipX " + (-l) n e iqX ". 

For the derivation of this expansion see [16]. Note that constants K n are exactly same as 
the corresponding ones in [10, 16]. 

Now all the statements in Theorem 1.2 follow from (2.1) and another key point: 3n > 2 
ct n K n 7^ 0. This is due to a contradiction argument. Assume that the assertion is not 
true, that is, Vn > 2 a n K n = 0. 

Case A. First we consider the case when B/ir is irrational. It is easy to see that K n ^ 
for all n > 2. Thus, a n = and this yields u(r,9) = aiJo(kr) for < r < rj and 
< 6 < ©. Since this right-hand side is an entire solution of the Helmholtz equation, 
the unique continuation property of the solution of the Helmholtz equation yields u(x) = 
cci Jo(k\x — x \) in R 2 \ D. This implies that u has to be bounded in a neighbourhood of 
y. However, since 

^o(x,y) ^ —log- 



2n \x — y\ 

as x — > y and E£>(-,y) is smooth in a neighbourhood of y, one knows that u(x) = 
&o(x,y) + Erj(x,y) is not bounded in any neighbouhood of y. Contradiction. 
Case B. Next consider the case when O/V is rational. One can write 

1 b 

- = 1 + 

7r a 

where a(> 2) and 6(> 1) are integers and mutually prime. Then we have 

{n > 2 | K n = 0} = {1 + l{a + b) \ I = 1, 2, • • •}. (2.2) 

Note also that \i+i( a +b) — al, I — 1,2, From the assumption of the contradiction 
argument one knows if n satisfies K n ^ 0, then a n = 0. From this together with (2.2) we 
have 

oo 

u(r, 9) = J2 &i+i{a+b)Jai{kr) cos aW. (2.3) 

1=0 

Hereafter we take two courses corresponding to Theorems 1.2 and 1.3. 
2.1 Completion of the proof of Theorem 1.2 

Since each al in (2.3) is an integer, the right-hand side of (2.3) gives a continuation of u 
onto (R 2 \(DU {y}))U B v (xq) as a solution of the Helmholtz equation and the continuation 
which we denote by u satisfies the rotation invariance in B rj (xo): 

u(r,0 + ^j=u(r,6). (2.4) 

Now having (1.1) and (2.4), one can apply Friedman-Isakov's extension argument [6] 
to u. See also [18] for the detail of the argument applied to a penetrable obstacle case. 
As a result one gets a continuation of u onto R 2 \ {y} as a solution of the Helmholtz 
equation. Since u(x) = u(x) = & (x,y) + E D (x,y) in R 2 \ D and En{x,y) is smooth in 
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a neighbourhood of y, one concludes that Ed( ■ , y) can be continued as a solution of the 
Helmholtz equation in R 2 . The continuation satisfies the Sommerfeld radiation condition 
and therefore has to be identically zero. This gives, in particular, u(x) = § (x,y) in 
B^Xq) and from (2.4) one gets 

$0(2:0 + rz(9),y) = $ (x + rz (o + -^j , yj , < r < rj, 9eR, (2.5) 

where z(9) = cos 9 a + sin 9 a L . 

Since both sides of (2.5) satisfy the same Helmholtz equation in \x — x \ < \y — x \, 
the unique continuation property of the solution of the Helmholtz equation yields: (2.5) 
is valid for all r with r < \y — Xq\. Choose a 9q in such a way that y = Xq + \y — Xq \ z(9 ). 
Since 2n /a < n, we have y ^ Xq + \y — x \ z(9q + 2ix / a). Then letting 9 = 9 and r ^ \y — x \ 
in (2.5), we have a contradiction because of the singularity of &o(x, y) as x — > y. 

This completes the proof of Theorem 1.2. 

□ 

Remark 2.1. In the proof of Theorem 1.1, we never make use of (2.4) after having (2.3) 
and instead take another course for the total field u = u( ■ ; d, k) in Theorem 1.1. 
The argument is as follows. Set 

di = cos 9 a + sin 9 a ± \o = Q- 7T , $i = cos 9 a + sin 9 a ± \ e= ^_ 7r ) +7r / 2 , 

(2.6) 

d 2 = cos 9a + sin 9a ± \e=- K , $2 = cos 9 a + sin 9 a ± | 6 ) =7r+7r / 2 . 

Note that d\ and d 2 are directed along the two sides that meet at xq. 

From the right-hand side of (2.3) one gets: for all r with < r << 1 Vu(io+rrfi)-t?i = 
and Vu(xq + rd 2 ) • $2 — 0. Then a reflection argument in [1] yields that this is true for 
all r > 0. However, from this together with the asymptotic behaviour of Vw ~ \/ e lkx - d as 
r — > 00 one gets d ■ $1 = d ■ <d 2 = 0. Contradiction. 

The advantage of this argument is: one does not need to use (1.1). In the following 
subsection we employ this argument after (2.3). 

2.2 Completion of the proof of Theorem 1.3 

We use the same notation as (2.6). First we claim that, as r — > 00, 

V^oix^y) ■ = l ^^e™/ 4 e ik ^-y> d ^u(y; -d 3 , k,x ) + O (J^ , (2.7) 

where Xj = x Q + rdj. 

This is proved as follows. Total field & D (- ,y) has the expression 



$ D (x, y) = %(x, y) + J qd ^y$o02, x)$d(z, y)dS(z), xeR 2 \D. (2.8) 

Since 



VMx,y) = ^^-XH^)'(k\x-y\), 
4 \x y\ 
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we have 

V,<I>„U, y) ■ 0, = ^^L^i {H Wy {k \ Xj _ y \). 

Here we note that (xj — y) ■ dj = (x — y) • dj + rdj ■ dj = (x — y) • dj since dj ■ dj = 0. 
This gives 

V x %(x 3 , y) ■ dj = J ^Jy^ iH^iklx, - y\). 
By (4.03) on p. 238 in [27], we know that as r — > oo, H^\r) and its derivatives satisfy 
if«(r) = ^| e -W4 r -i/2 e *r ^ + 0(r -i } ) 

(H^)"(r) = i^ e ™ /A r- 1/2 e ir (l + O^- 1 )) . (2.10) 
Using 

kj - y\ = r + (x - y) ■ dj + 0(r _1 ), 

and (2.9), we have 

(H^)'(k\xj - y\) = ^| e iT/4j b -i/2 r -i/2 e ifcr e ifc( a! o-y)-d i (! + ( r -i 



and 



and thus 



V^ fe,y) • ^ = ^ J/] ' ^ Jl^/^ e ik {xo -y). d] + ( r -5/ 2) _ (2.H) 



Let z G <9-D. We have 



$ (z,x) = -- v ■ ■ (#o 7(*k-*l) 



<9z/(z) ' 4 |x — 2 

and thus 



\dv(z) ' J A \x — z\ 
it (*-,). 1/(2) (x - ,) ^ _ ^ (»-^).^)(»-^) w _ 



4 |x — z\ a 4 |rr — 2 

This together with (x — z) ■ dj = (x — z) ■ dj yields 

d ^ f \ 1 q ik v{z) • dj 



Mz, \x= Xj ■ = -t i (C)'(^j - A 



^ 2 - Z) ■ u(z) (X -zl^( H (l)y,( k \ x . _ 



\Xj z\ 
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Note that the second term of this right-hand side is estimated by 0(r 5 / 2 ). It follows 
from these and (2.10) that 



v(z) ■ $j + ikdj ■ v{z) (xo — z) • dj / ^ ^/^^^{xo-zydj _|_ Q( r - 5 / 2 \ 
Now from (2.7), (2.8), (2.11) and (2.12), we obtain 

„-t,l/2 fa ikr i 

where 

U(y;d,k,x ) = (xo-y)- V M " K 
- / • 7?j + fifed, • 1/(2) (a* - z) ■ e ifc ( x °- z >^$ D (z, y)dS(z). 

Define 
Since 
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(2.12) 



(2.13) 



(9z/(z) 

one can rewrite (2.13) as 



((ar - ^) • #j e - ***"*') = - • t?j + ikdj ■ v(z) (x - z) ■ fy) e" 



■ikz-djU 



U(y;d,k,x ) = *j(y) + Q^fc) ' ^D^y)dS(z). (2.14) 

On the other hand, a combination of Green's identity, the Sommerfeld radiation condition 
for Wj(z) = w(z; —dj, k,xo) and &E>(z,y) and the boundary condition in (1.2) gives 



Wj 



dv(z) 



= I ^ r r{(xo-z)^ j e- ikzd ^ D (z,y)dS(z). 

JdD OV\Z ) 

Therefore we see that the left-hand side of (2.14) coincides with u(y; —dj, k,x ). This 
completes the proof of (2.7). 

Now the proof of Theorem 1.3 starts with having (2.3). By the same reason described 
in Remark 2.1, from (2.3) we have 

V x $ D (x j ,y)-& j = 0, (2.15) 
10 



where Xj = xq + rdj and < r << 1. First consider the case when y ^ xq + rdj for all 
r > 0. In this case a reflection argument in [1] ensures that (2.15) is valid for all r > 
and thus (2.7) yields 

u(y; —dj, k, x ) = 0. (2-16) 

If y = Xo + \y — Xo\dj, then V x §o(x, y) • $ j — for x = x + rdj with < r < |y — rr | 
and r > \y — x \. Then form (2.15) we have V x En{x,y) • $ j — for x = xo + rdj with 
< r << 1 and this is true for all r > by a reflection argument in [1]. Therefore we 
again have (2.15) for all r > \y — x \ and thus (2.16) too. 

Summing up, in any case we obtain equation (2.16) for j = 1,2. This is a contradiction. 

This completes the proof of Theorem 1.3. 

□ 

Remark 2.2. However, (2.16) is coming from only the leading term of the asymptotic 
expansion (2.7). Thus our next problems in this direction are as follows: 

(i) determine the complete asymptotic expansion of e~ lkr V x ^o{.Xj, y) ■'dj as r — > oo: 

oo 

e-^V^D^y)-^ ~ £ A m (y,£o,^K (3/2+m) - 

m=0 

(ii) If A m (y, x , dj, k) = for all m, then what happens? 

The main obstruction in this approach is the complexity of computing the asymptotic 
expansion in (i) as can be seen in the proof of Theorem 1.3. 



3 Other applications 

In this last section, instead we give three applications of the argument done in Theorems 
1.2 and 1.3. 



3.1 Thin obstacle 

It should be pointed out that the advantage of the assumption in Theorem 1.3 is that the 
result can be extended to a thin obstacle case. 

First we review a result in [16] which employs a single plane wave as an incident wave 
and corresponds to Theorem 1.1. 

Let £ be the union of finitely many disjoint closed piecewise linear segments denoted 
by Si, E 2 , • • ■ , E m . Assume that there exists a simply connected open set D such that D 
is a polygon and each Ej consists of sides of D. 

We assume that D C B R with a R > 0. We denote by v the unit outward normal on 
dD relative to B R \ D and set z/ + = v and v~ = — v on E. Given k > and d £ S 1 let 
u = u(x),x G R 2 \ E, be the solution of the scattering problem 

(A + k 2 )u = in R 2 \ E, 



= onE, 



lim Jr [ — ikw) = 0, 

r— >-00 J 
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where w = u — e tkx ' d , u + = w| R 2\j5 and u = u\ D . Note that this is a brief description of 
the problem and for exact one see [16]. Define 



I^(r;u,d,k) = \q^ Vt ~ ~q^ u ) dS ' 



In [16] we have established the following result. 

Theorem 3.1([16]). Letoo be regular with respect to E. If every end points of Hi, E 2 , • • • , E r 

satisfies x ■ oj < h^ioj), then the formula 

lim - log \Iz(r; u, d, k)\ = hz(u), 

T >00 f 

is valid. Moreover, we have the following: 

if t > h^(u), then lim T >oc e~ Tt \lY,{r; u, d, k)\ = 0; 

if t < hs(uj), then lim T >OQ e~ r *|is(r; u, d,k) \ = 00. 

If there is an end point x of some Ej such that x ■ ui = h-z(ou), then, for d that is not 
perpendicular to v on Ej near the point, the same conclusions as above are valid. 

Note that v on Ej fl B v (xq) for sufficiently small 77 > becomes a constant vector if 
Xq is an end point of E.,-. 

Here we present a result in which, instead of a single plane wave we make use of a 
single point source as an incident wave. 

Let y G R 2 \ D. Let E = EY,(x,y) be the unique solution of the scattering problem: 

(A + k 2 )E = inR 2 \E, 
3 -E^ = -^-%(-,y)onH, 



dv ± dv ± 

(BE \ 
— ikE\ =0. 
or J 

The total wave outside E exerted by the point source located at y is given by the 
formula: 

$s(z,y) = $o(x,y) + Ez(x,y), x G R 2 \ E. 

Given d £ S 1 choose i5 G 5 1 in such a way that -i? -1 = d. Let Xq G E and w = 
ws(i; — d, k,xo) be the unique solution of the scattering problem: 

(A + k 2 )w = inR 2 \E, 

{(x -x) ■#e- ihx d } on E, 



dv ± dv ± 

dr- 



iven. y/r(^- — ikw 



Define 

—d, k, x ) = (x — x) e~ lkx ' d + u>s(x; —d, k, x ). 
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Let R\ > R and y G dB Rl . Define 

J s (r; w, y, k) = jf ^ (J^®x(x, y) • v T (x; u) - -^ V A X ', w ) • y)^j dS. 

The following theorem is what we call an extension of Theorem 1.3 to thin obstacles. 

Theorem 3.2. Let u be regular with respect to E and let x G £ be the point with 
xo ■ oj = h-s(uj). Assume that y G satisfies uj;(y; —d, k, xo) ^ for a direction d E S 1 

that meets at x along a Ej . T/ien the formula 

^hr^-log |J s (r;u;,y, fc)| = 

zs vaZzc?. Moreover, we have: 

if t > hv(uj), then lim r >oc e~ r '| Je(t; y, k) \ = 0; 

i/t < hj:(uj), then lim T >OQ e~ rt \ Je(t; w, y, fc)| = oo. 

The proof of Theorem 3.2 is based on the convergent series expansion of $s( • , y) at a 
corner or end point of E. See Propositions 4.4 and 4.5 in [16]. Since the proof of Theorem 
3.2 can be done along the same line with that of Theorem 1.3, we omit the description. 



3.2 Obstacle in a layered medium 

We consider a medium that consists of two parts. One is given by R 2 \ Br and another 
is B R . We assume that the propagation speeds of wave in two parts can be different from 
each other. An obstacle D is embedded in B R as before. We assume that D is polygonal. 
Let us describe a mathematical formulation of the problem. 
Define 

' 7+, x G R 2 \ B R , 

7(x) = < 

k 7_, x G B R , 

where 7± are known positive constants. 

Fix y G R 2 \ Set £;+ = fc/ > /7+- Define 

®-r(x, y) = $+(x, j/) + e 7 (x, y), x G R 2 

where 



and e = e 7 solves 



$ + (x,y) = i ^(fe 4 .|x-y|) 



(V • 7 V + k 2 )e = -V • (7 - 7+) V$+(x, y) in R 2 , 

lim y/r ( — ik + e\ =0. 

r — >oo yy r j 
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Let E = E D:J (x, y), x G R 2 \ D solve 

(V • 7V + k 2 )E = in R 2 \ U, 

dE d 
7— = -7— $i(x,y) ondD, 

(dE \ 
- ik + E) = 0. 
or J 



Define 

$d, 7 (x, ?/) = $ 7 (x, y) + E D:J (x, y). 

Note that the existence and uniqueness of the solutions e and E can be established by 
using a variational formulation, for example, see [8]. 
Define 

K(r;u,y,k) = J (7. Jj-S^x, y) -v~{x;u) -^-^v~(x;u) ■ ^(x,y)\ dS(x), 



(3.1) 



where 



v;(x;oo) = e x<TUJ+i V^ TuJ± \ 
$B, T (z, V) = $d, 7 (z, y), xeB R \D 

and A;_ = kj ^/T^- 

Theorem 3.3. Assume that oj is regular with respect to D and that 

diamD < dist(D,dB R ). (3.2) 

Moreover assume that there exists a j e {l,---,m} snc/i that A; 2 noi a Neumann 
eigenvalue for — A m Dj. It holds that 

lim^^- log \K(r-,u,y,k)\ = h D (u). 

Moreover, we have the following: 

if t > hoioj) , then hm e~ Tt \K(t; oj, y, k)\ = 0; 

if t < hoioj) , then lim e~ Tt \K{r; oj, y, k)\ = 00. 

Proof. Instead of u(x) = $d(x, y) in the prof of Theorem 1.2 set u(x) = <&£,^(x,y). For 
this u the same argument with (3.2) instead of (1.1) as described in Cases A and B in the 
proof of Theorem 1.2 works. Thus we have a continuation u of u in B R \ D onto B R as a 
solution of the Helmholtz equation Au + k 2 _u = in B R . Since u satisfies the Neumann 
boundary condition du/dv = on dD, u also satisfies the condition on 3D and by the 
assumption k 2 _ , it must hold that u = in Dj for some j. Then, the unique continuation 
theorem for the Helmholtz equation yields u = in B R \ D. This yields that the Cauchy 
data of $d i7 (ie, y) on dB R vanish and thus $D i7 (x,y) = for x G R 2 \ (B R U {y}) by 
the uniqueness of the Cauchy problem for the Helmholtz equation with wave number k + . 
Since $d i7 (x, y) is singular as x — > y, this is a contradiction. 
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□ 

In this theorem the data are given by the Cauchy data of the total wave field on dB R . 
It is an interesting open problem when the receivers are located on dB R+e with a e > 
how one can apply the enclosure method in an explicit form. 

Here we propose one heuristic approach based on Theorem 3.3 in the case when e is 
sufficiently small. 

Assume that we have <&d j7 (x, y) for all x G dB R+e exactly. Solve the exterior problem 
in R 2 \ B R+e : 

+ k 2 + ^> = inR 2 \B R+e , 
V(x) = $ D „(x,y) - $ + (x,y) ondB R+e , 

lim y/r — ik + ^ = 0. 

r — >oo yg r j 

Then we have V(x) = ® Djl (x,y) - & + (x,y) for x G R 2 \ B R+e . This gives 

We use for the computation of the Cauchy data of $d i7 (x, y) on <9-B# from outside B R 
the approximation: 

®d,i(x, y) « $D, 7 (a; + ei/(x), y), 

<9$ D7/ . 9$+, , , , 9*, 
g ' {x, y) ~ -^-(^ + ev \ x )iy) + ~frj( x + ev \ x ))- 

Using these computed Cauchy data from outside B R and the transmission condition 

®D,~f(x,y) = $D !7 (x,y), 7+^^(2;, y) = 7_^?^(x,y), x G dB R 

which is implicitly included in the governing equation, we compute K(r;u,y,k) by re- 
placing $^ 7 (a;, y) and 7_(<9$£ /du)(x, y) in the right-hand side of (3.1) with <&D n (x + 
ev{x) 1 y) aiid ^ + {(d$ + / dv)(x + ev(x) ,y) + (d^! / du)(x + £u(x))} , respectively. Clearly, the 
effective range of r shall depend on the size of e. 

It would be interesting to test this approach numerically and check its performance. 
This belongs to a next research plan. 

3.3 Unsolvability of the far-field equation for polygonal obsta- 
cles 

Let k > and d G S 1 . Let F D ((p; d, k) denote the far- field pattern of the scattered wave 
w{x) = u(x; d, k) — e lkx ' d . 

Given y G R 2 the far-field equation for unknown g G L 2 (S l ) 

/ F D ( V - d, k)g(d)dS(d) = —= e -^y, ip E S 1 (3.3) 
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plays the central role in the linear sampling method [3]. 

Note that the right-hand side of (3.3) coincides with the far-field pattern of the field 
<&o(x,y) with x = rip as r — > oo; the left-hand side of (3.3) coincides with the far- 
field pattern of the scattered field w = w g which is the unique solution of the scattering 
problem: 

(A + k 2 )w = inR 2 \ D, 

dw dv 9 a n 
ov av 

lim Jr l^r- — ikw] = 0, 
where v g denotes the Herglotz wave function with density g: 

:) = J e ikx - d g{d)dS{d), x G R 2 . 



Vg(X) 



Note that w g satisfies w 9 \b r G H 1 (Br\D) for a sufficiently large R and the inhomogeneous 
Nuemann boundary condition on 3D should be considered in a weak sense. 

In this section, using the idea of the proof of Theorem 1.2, we give a proof of unsolv- 
ability of equation (3.3) for any k > and y G R 2 provided D is polygonal. 

Theorem 3.4. For any k > and y G R 2 there exists no solution g of equation (3.3). 

Proof. We employ a contradiction argument. Assume that equation (3.3) admits a 
solution g. Then the coincidence of both fa-field patterns of w g and <&o(-,x) yields 
w g (x) = $o(x,y) for x G R 2 \ Br with a sufficiently large R. From the unique continua- 
tion property for the Helmholtz equation this coincidence gives 

w g (x) = $ (x,y), Vx G (R 2 \D) \ {y}. (3.4) 

Since w 9 \ Br G ^(Br \ D) and $ ( • ,u')\b r does not belong to ^(Br \ D) for all y' G 
R 2 \ D from (3.4) one gets y G D. Note that this part or this type of argument is well 
known in the linear sampling method. It shows that if (3.3) is solvable, then y G D. The 
problem is the next to the intermediate conclusion y G D. Now we have 

w g (x) = $ (x, y), Vx G R 2 \ D. (3.5) 

Define 

Ug(x) = Vg(x) + Wg(x), X G R 2 \ D . 

Note that u = u g satisfies the Helmholtz equation in R 2 \ D and the homogeneous 
Neumann boundary condition dujdv = on dD. From (3.5) we have u g (x) = v g (x) + 
$o(x, y), x G R 2 \ -D and this right-hand side gives a continuation u g of u g onto R 2 \ {y} 
as a solution of the Helmholtz equation. 

Choose a cu G S 1 that is regular with respect to D and define 

'^=L& Vr -^ a °) dS ' T>0 ' m 

where v T (x) = e ^-(™+ i ^ 2 + fc ^ ± ). 
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Let xq € dD with xq ■ u = ho(oj). One may assume that y G D\, where D\ is a 
connected component of -D. Then one can choose a small 5 > such that if |x — y| < 5, 
then x & D 1 and rr-o; < h D (u) — 5. Replacing m 9 in (3.6) with u g and applying integration 
by parts, we obtain, as r — > oo 



(3.7) 



= e"^) / f - dS 1 = 0{re-^). 

J\x-y\=8 \ OV OV ) 

Hereafter we make use of the same notation as those of the proof of Theorem 1.2. 
Recalling boundary condition du g /du = on dD, one has the expansion 

oo 

u g (r,9) = ctiJv(kr) + a n J\ n (kr) cos \ n 9, < r < rj, < 9 < 
and applying the argument for deriving (2.1), we obtain 

/(r) e -»W,o^ e -r^H ~ -i ^ 6 fc An anJt " , (3.8) 

n=2 S 

where X n are constants given by the formula K n = e l P A - + (_i)"- e J <? A - anc i s — a/t 2 + P+t. 
Since (3.7) implies that e~ ThD ^I(T) is rapidly decreasing as r — > oo, all the coefficients 
of the right-hand side of (3.8) have to vanish, that is 

a n K n = 0, Vn > 2. (3.9) 

First consider the case when Q/tt is irrational. It is easy to see that fT n 7^ for 
all n > 2. Thus from (3.9) one gets a n = and this yields u g (r,9) = aiJo(kr) for 
< r < i] and < # < 0. Since this right-hand side is an entire solution of the 
Helmholtz equation, the unique continuation property of the solution of the Helmholtz 
equation yields u g (x) = a\ J (k\x — x \) in R 2 \ {y} and thus one gets 

%(x, y) = a^Jo^klx - x \) - v g (x), x^y. 

Comparing the behaviour as x — > y on both sides, we obtain a contradiction. 

Next consider the case when Q/tc is rational. Applying the same argument for the 
derivation of (2.4), we have a continuation of u g onto (R 2 \ D) U B^xq) as a solution of 
the Helmholtz equation and its continuation which we denote by u' satisfies the rotation 
invariance 

u (r, 9 + -^j = u'(r, 9), < r < r], 9 e R, 

where a > 2 is an integer. Since the unique continuation property gives u g (r, 9) = u'(r, 9) 
for < r < i] and thus one gets 



u g 



(r, 9 + ^p)= u g (r, 9), < r < r], 9 e R. (3.10) 
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Since u g satisfies the Helmholtz equation for \x — x \ < \x — y\, it follows from the unique 
continuation property and the rotation invariance of the Helmholtz equation that r] in 
(3.10) can be replaced with \x — y\: 

u g (x + rz(9)) = u g (x + rz(e + , < r < \x - y\, 9 G R, (3.11) 

where z{9) = cos 9 a + sin 9 a L . Now choose a 9q in such a way that y = x + \y — xq\ z{9q). 
Since 2n/a < n, we have y ^ Xq + \y — x \ z(9 + 2n/a). Then letting 9 = 9 and r f |y — x \ 
in (3.11), we have a contradiction since = v g (x) + $ (£, y) for a; 7^ y and 

1 1 



2/T |x — y| 

as x — > y. 
□ 

Using a variational formulation in, e.g., [8], one can formulate and establish the unique 
solvability of the scattering problem of acoustic wave by a sound-hard obstacle D with 
Lipschitz boundary. We use the same notation as those in the case when D is polygonal. 
Having the far-field pattern for D with Lipschitz boundary, one can extend Theorem 
3.4 to a slightly general case. We say that D with a Lipschitz boundary has a horn, 
if there exist a u E S 1 that is regular with respect to D and 5 > such that the set 
V = {xED\x-u> h,£>(oj) — 6} becomes a finite cone with the vertex at the point in 
{x I x ■ 00 = hr>(oj)} H dD and the base on x ■ u = hn{oo) — 5. We say that V is a horn. 

The conclusion of this section is the following statement and since the proof is really 
a minor modification of that of Theorem 3.4 we omit the description of the proof. 

Corollary 3.1. If D with a Lipschitz boundary has a horn, then for any k > and 

y G R 2 there exists no solution g of equation (3.3). 

Note that, in [26] the far- field equation for a single circular obstacle with an arbitrary 
radius has been considered and it is shown that the equation is not solvable except for its 
center point. Corollary 3.1 means that the existence of a horn V even it is small prevents 
the existence of solution of (3.3) for any k > and y G R 2 . 

It should be pointed out that the linear sampling method is not based on the solvability 
of the far-field equation. Instead a family of approximate solutions of the far-field equation 
is taken. See [2] for interesting study of the method itself. 
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